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In this paper, we propose opportunistic interference alignment (OIA) for three-transmitter multiple-input multiple- 
output (MIMO) interference channels (ICs). In the proposed OIA scheme, each transmitter has its own user group 
that consists of K users, and each transmitter opportunistically selects the user whose received interference signals 
are most aligned. Thus, three-transmitter MIMO IC is opportunistically constructed by three transmitters and their 
^2 \ selected users. Contrary to conventional I A, perfect channel information for all of the interference links is not 

0^ ' required at the transmitter Each user just needs to feed back one scalar value in the proposed OIA scheme. When 

o . 

the number of receive antennas is Nr (with Nr = 2M) and the number of transmit antennas is Nt (with Nt > M), 
we prove that each transmitter can achieve M degrees of freedom (DoF) as the number of users in each group goes 
' to infinity (K — > oo), in comparison with the DoF known to be achievable in the three-transmitter M x 2M 

MIMO IC. We also show that the user scaling K cx P"*^' enables the network to achieve M DoF per transmitter 



where P is transmit power, and a is the relative path loss of the interfering channel in decibels. The computational 
complexity of the proposed OIA scheme is analyzed and compared with conventional user selection schemes. 
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I. Introduction 

Interference alignment (lA) has been touted as a key technology for handling interference in future wireless 
communications. Contrary to the conventional wisdom that orthogonalization is the best way to handle interference, 
ifTl showed that lA can achieve a total of ^ degrees of freedom (DoF) in an A^-transmitter single-input single-output 
(SISO) interference channel (IC). Following this work, there have been many intensive studies on lA Hl-llTl. It has 
been recently shown that the achievable DoF at each transmitter in a A^-transmitter Nt x Nr MIMO IC becomes 

m 

DoF = mm(NT, Nr), N <R 

MIMO IC 



R + l MIMO IC R + 1 



ma.x{NT,Na) 
mm[NT,NR) 



. In 



where Nt and Nr are the numbers of transmit and receive antennas, respectively, and R 
the N > R case, symbol extensions are needed to achieve the DoF in time varying or frequency selective fading 
channels. Despite the promising aspects of lA, its implementation has many challenges. lA requires perfect global 
channel knowledge (including whole channel links) at the transmitter, which involves excessive signal overhead. 
Imperfect channel state information significantly degrades the gain of lA. The large computation complexity 
necessitated is also regarded as a big challenge for practical implementation. The sub-optimality of lA in the 
practical operating SNR region is another problem. 

Recently, lA techniques to ameliorate these problems have been investigated. Iterative lA algorithms were 
proposed to optimize precoding matrix and to reduce the global channel knowledge burden based on channel 
reciprocity 161. To reduce computational complexity, Suh and Tse Q proposed a subspace interference alignment 
technique for an uplink cellular network system. In ||8|, lA was opportunistically performed in MIMO cognitive 
radio networks, where secondary transmitters transmit their signals on only spatial dimensions not used by primary 
transmitters. lA with imperfect channel state information (CSI) was shown to achieve the same DoF as lA with 
perfect CSI if the feedback size per user is properly scaled lO, ifTOl . In ||9l, a DoF of ^ was shown to be achieved 
with limited feedback in an A^-transmitter frequency-selective SISO IC if the feedback size per user is maintained 
as N{L — l)logP bits, where L is the number of taps and P is the total available transmit power. lITOl showed 
that the maximum number of DoF is achievable with limited feedback in an A^-user MIMO IC if the feedback size 
scales with the SNR. Also, full DoF is known to be achievable with analog feedback if the analog feedback power 
grows with the transmit power lITll . 
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A. Opportunistic Interference Alignment ( OIA ) 

Although there have been significant efforts to overcome the practical challenges of lA, the inherent shortcomings 
of lA highly motivate the development of more practical lA techniques that work well with limited feedback and 
reduce computational complexity while maintaining the promised gain. In this paper, we propose an opportunistic 
interference alignment (OIA) scheme in a three-transmitter MIMO IC with a generalized antenna configuration. In 
our system model, there are three transmitters and three user groups corresponding to the transmitters. Only a single 
user in each user group is served by each transmitter, so a three-transmitter IC is opportunistically constructed. In 
our OIA scheme, the concept of opportunistic beamforming is incorporated in lA. Each transmitter transmits with 
random beams and each user feeds back the distance between the two received interfering channels. The feedback 
information corresponds to a measure of interference alignment and is relatively small. At each transmitter, the 
user whose interfering channels are most aligned is selected and served. The proposed OIA scheme in this paper 
generalizes our preliminary study on OIA |[T2l by adopting multiple streams at each transmitter. Compared to this 
paper, |[T2l just considered two antennas at each transmitter and receiver, respectively, and assumed only a single 
stream from each transmitter. 

Contrary to opportunistic beamforming in a MIMO broadcast channel lfT3l . lfT4l . random beams generated at a 
transmitter is not used for user selection at the transmitter, but helps other transmitters select their users. Once the 
users whose interfering channels are most aligned have been opportunistically selected, the selected users perform 
postprocessing with the equipped receive antennas. Because the postprocessing is performed after interference 
aligned user selection, only the selected user (as opposed to all users) needs to compute the postprocessing matrix. 
Correspondingly, the computational complexity is far reduced compared to conventional opportunistic user selection 
schemes where every user feeds back the signal-to-noise ratio (SNR), interference-to-noise ratio (INR), and the 
signal-to-interference -plus-noise ratio (SINR). 

B. Contributions 

We investigate the optimality of the OIA scheme in a three-transmitter Nt x Nr MIMO IC where K users are 
associated with each transmitter and the selected users together with their transmitters construct a three-transmitter 
IC. 

• If each receiver has Nr = 2M antennas and the transmitter has more than M antennas such that Nt > M, 
we prove that each transmitter achieves ^ = M DoF without symbol extension by the OIA scheme as the 
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number of users in each user group, K, goes to infinity. For a three-transmitter M x 2M MIMO IC, ^ DoF 
per user is known to be achievable (with symbol extension) according to the result given in (dJ. This result 
seems to to be contradictory at first glance, but the required spatial dimensions for M data streams are secured 
through the user dimensions provided by the K users. This means that multiuser DoF are translated into lA 
spatial dimensions. 

• We show that the number of users associated with each transmitter, K, should scale as oc P"*^^ to achieve 
M(= Nji/2) DoF per transmitter where a is the relative path loss of the interfering channel in decibels. 

• When K is finite, the achievable DoF of the OIA scheme is proved to be (1 — a)M. To achieve DoF M' such 
that (1 — a)M < M' < M, we show that the required number of users K should increase on the order of 

K o(. py " ) . 

• Finally, we look into the practical advantages of the proposed OIA scheme; it is shown that the proposed OIA 
scheme significantly reduces not only the amount of feedback information but also computational complexity 
while achieving a notable rate improvement compared to conventional opportunistic user selection schemes. 

C. Organization 

The rest of this paper is organized as follows. In Section II, our system model is described. Preliminaries about 
the angles between two subspaces are provided in Section III. The proposed OIA scheme is described in Section 
IV, and the achievable rate and DoF are analyzed in Section V. Several conventional opportunistic user selection 
schemes are summarized and compared with OIA scheme in Section VI. The conclusions and comments on areas 
of future interest are given in Section VII. 

D. Notations 

Throughout the paper, A^, Ara(A), v„(A), tr{A) and ||A||i? denote the conjugate transpose, nth largest eigen- 
value, eigenvector corresponding to A„(A), trace, and Frobenius norm of matrix A, respectively. Also, I„, diag{-), 
C" and C™^" indicate the n x n identity matrix, a diagonal matrix whose diagonal elements are (•), the n- 
dimensional complex space, and the set of m x n complex matrices, respectively. 

II. System Model 

As depicted in Fig. [T] there are three transmitters having Nt antennas, and each transmitter has its own user group 
consisting of K users with antennas each. In our system model, only a single user in each group is selected 
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and served by each transmitter. The transmitters and their selected users construct a three-transmitter MIMO IC. 
We assume that each transmitter sends M(< Nt) streams and each user has Nji{= 2M) antennas. The transmitters 
are assumed to have only partial CSI fed back from each user, and no information is shared among the transmitters. 
In our system, user selection and transmission take the following four steps: 

• Step 1: Each transmitter transmits with random beams. 

• Step 2: Each user feeds back one analog value. 

• Step 3: Each transmitter selects one user from its user group. 

• Step 4: Each transmitter serves the selected user. 

In the first step, the ith transmitter transmits its signal using a random beamforming matrix Wj S C^^xM ^^^j^ 
that w|Wj = Im- The received signal at the kth user in the ith user group is modeled as 

3 

3 

where H^^;^^ G (^NrxNt denotes the channel matrix from the jth transmitter to the A:th user in the ith user group 
whose element is an independent and identically distributed (i.i.d.) and circularly symmetric complex Gaussian 
random variable with zero mean and unit variance. The transmit signal from the ith transmitter is denoted by 
Xj G C^^^^. Because a transmitter does not have CSI for power allocation, the same power is allocated for each 
stream and the transmitter power constraint is given by E{xjxJ} = ;^Ia/- Hj ^ G (j^N^xi denotes Gaussian noise 
following a complex Gaussian distribution with zero mean and an identity covariance matrix rij ^ ~ CA/'(0, Iat^ ). 7 
is a relative path loss from the interfering transmitter to each user and is assumed to be the same for all transmitters 
and users. We denote the received average interference power from each interfering transmitter as Pj and define 
a constant a which is the relative propagation path loss of the interfering channel in decibels. Then, the received 
interfering power Pj becomes Pj = 7P = P", and it is satisfied that 7 = P"~^. We also assume the average 
interfering power is less than the desired signal power such that Pj < P (equivalently, < a, 7 < 1). The effective 
channel from the jth transmitter to the kth user in the ith user group becomes H^-''^ G C^^^^'^. Because a random 
beamforming matrix Wj is unitary and independent of H >^ , each element of H^^ ^ becomes an independent and 
identically distributed (i.i.d.) complex Gaussian random variable with zero mean and unit variance. At each user, 
the received signal can be equalized using the postprocessing matrix. We assume the postprocessing matrix at the 
kth user in the ith user group, denoted as Fj ^ G C^^^^^, satisfies Fj ^Ft^ = Then, the received signal after 



postprocessing at the user becomes 



and the achievable rate of the user is given by 



log2 



log2 



P 



M 



(i)TT(i)tT^t 



Im + Fj^fc ( ^H^*fcH^*fc + EjVi j F[^ 



(2) 



In Step 2, each user feeds one analog value back to its own transmitter. The type of information needed to be fed 
back for user selection varies with the type of postprocessing and user selection scheme. In the OIA scheme, the 
chordal distance between the two subspaces spanned by the two interfering channels is used as feedback information. 
The details of feedback information construction will be given in Section IV. In Step 3, each transmitter selects a 
single user among the users in its user group. After user selection at each transmitter, the selected user is served by 
the random beams. If the kith user is selected in the ith transmitter, the average achievable sum rate of the system 
is obtained by 

3 
i=l 

III. Preliminaries - Angles between Two Subspaces Spanned by Two Interfering Channels 

The Stiefel manifold and Grassmann manifold are widely used geometric concepts in communication systems 
HS-IIl. The Stiefel manifold SNn,M{Q (for Nr > M) is the set of all unitary matrices in C^"''^^ i.e.. 

If the columns of two matrices Si, S2 G SNa,M{C) span the same subspace. Si and S2 are defined to be equivalent 
lITSl . As the quotient space of 5iv^,j\/(C), the Grassmann manifold QNa,M{C) is the set of all M-dimensional 
subspaces in an A'^/j-dimensional space C^". 

Consider two M-dimensional subspaces A,B in C^", i.e., A,B€^ Qnr,m- define the angles between the 
subspaces, we use the principal angles, also called the canonical angles. The principal angles 9i, . . . , 6m G [0, 7r/2] 
between A and B are recursively defined by searching for A'^/j-dimensional vectors {aj} and {bj} such that |[20l 
Chap. 12] 



cos^,- = max max la^bl = lajb,- 
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subject to ||a|| = 1, ||b|| = 1, a^'a^ = 0, b^bj = (1 < j < i — 1). The vectors {aj} and {bj} corresponding to 
the principal angles are called principal vectors of A and B, respectively. 

The distance between subspaces can be defined in many ways and chordal distance is arguably the most widely 
used one among them. The squared chordal distance between A and B (denoted by d1{A^B)) is represented with 
the principal angles by 

M 

dl{A,B) = ^s\n^ Oi. (3) 

Alternatively, the squared chordal distance can be represented with generator matrices; A G ([^NrxM g ^ 
(^NbxM ^j.g called generator matrices if their columns are orthonormal, respectively, i.e. A^A = B^B = Im, 
and span the subspaces A and B, respectively. Thus, the generator matrices belong to the Stiefel manifold, i.e., 
A, B G 5Ar^^A/(C). Although generator matrices A and B are not unique, the generator matrices for A and B are 
equivalent in the Stiefel manifold, respectively. Also, the squared chordal distance between A and B is uniquely 
defined with arbitrary generator matrices A and B such that 

dl{A,B) = ^\\AA^-BB^\\l 

= M-tr(Al'BB^A). (4) 

The principal angles and principal vectors can also be obtained from the generator matrices. Let the singular 
value decomposition (SVD) of A^B be denoted by ll20l Chap. 12] 

A^^B = YDZT (5) 

where Y, Z G ([^MxM unitary matrices and D = diag{^i, ^2, ■ ■ ■ , fJ-Ai) such that fii > fi2 ■ ■ ■ > fJ'M > 0. Then, 
the ith singular value, fii, is given by 

fii = cos 6i 

where 9i is the ith principal angle and corresponding principal vectors a^ and bj become, respectively, 

aj = Ayi, bj = Bzj 

where and Zj are the ith column vectors of Y and Z, respectively. 

In our system model, each user has two interfering channels, and each interfering channel spans M-dimensional 
subspace in C^^ with probability 1 because the elements of each A'^R-dimensional interfering channel vector are 



i.i.d. and circularly symmetric complex Gaussian random variables. Thus, the union of two subspaces formed by 
two interfering channels is a 2M{= A'^r) -dimensional space with probability 1 for the same reason. 

Lemma 1: If two subspaces A and B are subspaces formed by interfering channels, the eigenvalues of AA^^+BB^ 
can be represented with the principal angles such as 



V ' V ' 

M M 

where A and B are arbitrary generator matrices of A and B, respectively. 



(6) 



Proof: See Appendix lAl 



IV. Opportunistic Interference Alignment 

This section describes the proposed OIA scheme. Without loss of generality, we consider the first transmitter 
and the users in the first user group. From Q, the achievable rate at the kth user in the first user group is given by 



Cifc = log 



We decompose ([7]) into the rate gain term and the rate loss term denoted as C^^^" and C}°^^ respectively, given as 

i„+F,4iiHs;x? + ^EHS:lH«')Ft, 



(7) 



^gain 1 
CU =log2 



i=2 



pa 



M 



i=2 



(8) 
(9) 



and it is satisfied that Cik = C- 



gain 
l,fc 



A. Feedback Information from Each User 

If the interfering signals are approximately aligned in the M dimensional space, the residual interference in the 
desired signal space will be negligible. The level of interference alignment can be effectively measured by the rate 
loss term after postprocessing when the user with the minimum rate loss is selected in the OIA scheme. If we 
assume that each user finds the postprocessing matrix to minimize the rate loss term, the minimum rate loss term 
at the /cth user in the first user group becomes 



mill C^°j^ = mill log2 



pa 



i=2 



and the minimum occurs when 



Fl,A: = [vj\/+l(Bi^fe), . . . , V2Af (Bi^fc)]^ 



(10) 



where Bi ^ = Y^^=2^i\^i^k ■ '■^^ ^^^^ minimum rate loss term is selected at the transmitter, the 

rate loss term at the selected user becomes 



min min Cff = min min loEfa 
k F' ^•'^ k F' 



pa 

IvI 



i=2 



In average sense, this rate loss term at the selected user is upper bounded as 



■jj(2) jj(3) 



mm mm Oi u 
k F' 



(a) 

< EvV(2) VV(3) 

"■l,kJ"-l,k 

(b) 

< E£j(2) -£^(3) 



min min E 



P 



2M 



mm Mlog2 ( 1 + ^ XI -^"^ i^hk) 



m=M+l 



(i) (i) ~ (i) 

where A\ ^ is a diagonal matrix comprising unordered eigenvalues of fjl\ ^ , ^ is a unitary matrix constructed 



with corresponding eigenvectors, i.e. Hj^ ^H^ ^ 



and Ci,, = T.U^%^?1 In the above 
equations, inequality (a) holds because h[*|, and A^i \ are independent [1191, and the average of the minimum 
values is smaller than the minimum of average values. Also, inequality {h) is satisfied because 



min E.( 
F' 



rnin E.(2) .(3) logo 

'^l,k>'^l,k 

pa 



pa 



1=2 



< min logn 



(d) 



< 



lM + i^F'(5;H«E[A«]H«t)F't 

i=2 

logs n (l + ^"A^(Ci,fc)) 

m=A/+l ^ ^ 

/ pa \ 
Mlog2(l + — A^(Ci,fc) 

^ m=A/+l ^ 



(11) 

(12) 
(13) 



where E[a[*|j] = MIm l[T9l . In above equations, inequality (c) is from Jensen's inequality, and equality (d) is 
when Fi = [va/+i(Ci fe), . . . , V2M(Ci^fc)]^ Also, inequality (e) is from the concavity of logarithm function using 
Jensen's inequality. 

In order to avoid the large burden required by having all users compute their postprocessing matrices given by (fTOl) 
for feedback and user selection, each user feeds back only the upper bound on its minimum rate loss given in (fT3l) 
instead of the minimum rate loss itself. Correspondingly, each transmitter selects the user with the smallest rate loss 
bound. Note that the rate loss bound given in (fT3l) does not depend on postprocessing allowing the computational 
complexity can be significantly reduced at each user. Consequently, the required feedback information from user k 



associated with the first transmitter is given by 

2AI 



\ /'H(2)VT(2)t , fT(3)TT(3)t 



(14) 



m=M+l 
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and the following theorem shows that the required feedback information is equivalent to the chordal distance 
between the two interfering channels. 

Theorem 1: The required feedback information for the OIA scheme given in ([T4l ) is equivalent to the chordal 
distance between the two subspaces spanned by the two interfering channels such that 

2M 
m=M+l 



Proof: From Lemma \T\ the eigenvalues of H^^],h[^],^ + h[^].h['^|,^ in descending order become 

1 + li, . . . , 1 + tiIj, 1 - ijIj, ...,l-r]l, 

^ V ' V ' 

M M 

where ?7i > . . . > tjm correspond to the cosine of principal angles (/>!,..., 0a/ between the two subspaces generated 
by the columns of H.f'l and h[^], (i.e., r/j = cosc^j). Thus, the feedback information given in (fT4b becomes 

M M 



fi /u{2) £r(3) 

1=1 i=l 



from the definition of chordal distance given in (jSj. ■ 
Thus, from the definition of chordal distance given in the required feedback information from the fcth user 



in the first user group becomes 



<'?(H!lHS) = M-ir(H™HS>S'HS5: 



where h[^].^ and H^"^].^ are the generator matrices of h[^]. and H^"^]., respectively, and can be found via SVD or 
QR decomposition. 



B. User Selection at the Transmitter 

The first transmitter selects the user k!^^^ who has the smallest feedback value such that 

fc™^ = argmin dl{Ylfl,,Ylfl). (15) 

l<k'<K 

Define a random variable D as the feedback value of the selected user which is equivalent to the chordal distance 
between the interfering channels of the selected user, i.e., 

D = dc(H^^^.oiA, H^'^^oia)- (16) 
Then, the average value of D is given in the following lemma. 
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Lemma 2: The average feedback value of the selected user, E[Z)], is equivalent to the quantization error to 
quantize H G £_NrxM averaged over all possible random codebooks of size K such that 



mO] =E min (f^(H,W) 



(17) 



where Crand C ^iv^,j\/(C) is a random codebook of size K whose codewords are independently generated from the 
isotropic distribution on Qnr,m{^)- 

Proof: See Appendix |B] ■ 

It has previously been shown that the average quantization error to quantize a randomly distributed source on 
the Grassmann manifold QNuMi^) ^^^^ the random codebook Ciand C t/Ar„,A/(C) of size K is upper bounded by 
D such that Wl 



E 



min d (H,W) 

WeCrand 



where D is given by 



M2 



(18) 



with 7] = p^j/a^.]^) rii^i ^r{M-i+i) ' ^ ^ ^ ^^^^ number chosen to satisfy {r]K)^ < 1. Thus, we can 

conclude that the average feedback value from the selected user is upper bounded as 



E[D] < D. 



(19) 



Note that the second term in ([TSl l can be negligible for large K ifTTl . For most engineering cases, the main order 
term of (fTSl ) is sufficiently accurate |[T7l-|[T9l. 



C. Postprocessing at the Selected User 

As noted in Theorem [T] each user computes only the chordal distance between two interfering channels and 
feeds it back to its transmitter. Only the selected user, A;p'^, applies a postprocessing matrix to minimize the rate 
loss term in (fTOl) . The postprocessing matrix at the selected user becomes 



- VA/+l(BfcOiA), • • • ,V2m(B;,oia) 



(20) 



f(2) „(2)t , tt(3) TT{3)t 



It should be noted again that the complexity of the proposed OIA scheme can be significantly reduced because 
only the selected user applies a postprocessing matrix. The complexity of the OIA scheme is analyzed in the latter 
part. 
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V. Achievable Rate and Degrees of Freedom (DoF) 

This section analyzes the achievable rate of the proposed OIA scheme and its DoF. Without loss of generality, 
the average achievable rate and a DoF of the first transmitter are derived. To simplify notations, we omit the user 
indices in the variables such as 



where Fi is the postprocessing matrix of the selected user given in (1201 ). Using the simplified notations, the 
received signal after postprocessing at the selected user at the first transmitter is rewritten by 

3 

Fy = FH^^^xi + FH^^^Xi + Fn. 

i=2 

When each group has K users, we denote the average achievable rate of the selected user using OIA scheme 
as R^^, and the average rate gain term and the average rate loss term at this time as -R^™ and -R|°|], respectively. 



i.e., R^^^ = fco,A, i?"^" = EC^'^oiA, and i?}'^^ = ECl^^U- Using the simpUfied notations, R'^^, Rf^^ and i?]^^^ 
are represented as 



Epi-f = Elog2 



i?fr;? = Eiog2 



Im + ^FH(i)H(i)tFt + ^ FH«H«tFt 



=Elog2 



lM + ^ELFHWH«tFt 



1=2 



i=2 



(21) 
(22) 



so it is satisfied that = R^^^ - R^^y 



Theorem 2: When the number of users, K, is finite, the achievable DoF of the proposed OIA scheme becomes 
(1 - a)M such that 



lim 



pOIA 



(1 - a)M. 



(23) 



P^oo log2 P 

Proof: We can directly derive the achievable DoF from (|2TI ) and (l22l ). When K is finite, it is satisfied that 



lim 



P— >-oo 



M because the matrix scaled by P is dominant compared to the matrix scaled by P° (0 < a < 1). 



IK] 



Also, it is satisfied that limp_>oo p 



aM because H^^^H^^^^ + H^'^^H^'^)^ has non-zero M eigenvalues under 



the finite number of users with probability one. Thus, ( |23] | is hold from R'^^ = ^^^"i ~ ^\ky ' 

Fig. [2] shows the average achievable rate of each user with the proposed OIA scheme according to a when 
K = 10. As stated in Theorem [2l the achievable DoF of each user becomes (1 — a)M when the number of users 
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K is finite. When a = 0, the achievable rate linearly increases with M so that full DoF M is achieved only with 
finite K. When a = 1, the achievable rate is saturated in the high SNR region and hence the achievable DoF 
becomes zero. 

Lemma 3: When the number of users in the each group is K, the average rate loss term of the selected user, 



is bounded as 



where D is given in (1181) . 

Proof: See Appendix O ■ 
Theorem 3: When the number of users per group goes to infinity [K — > oo), the achievable rate of the selected 



user, R?^^, behaves like 



|oo, 

1 



(2P) 

where C^^^ „^ xAi^'^^ P/M) is the achievable rate of the block diagonalization system where a transmitter with 
Nn antennas serves ^ users with M receive antennas each, and M independent streams are transmitted for each 
user with power per stream under perfect CSIT. 

Proof: See Appendix iDl ■ 

Corollary 1: The achievable DoF of the selected user becomes M as the number of users goes to infinity such 
that 

^OIA 

lim lim — = M. (24) 

P-s-oo K-^oo log2 P 

Proof: Since the achievable DoF of a BD system is given by 

lim ^2X2xA/(2^) ^ 

P--s>oo log P 

the equation (l24l ) is easily obtained from Theorem [3] ■ 
Because the achievable DoF of each user in conventional lA is known to be at most ^ when Nr = 2Af and 
Nt = M, Theorem |2] and Corollary [T] appear contradictory at first glance. In case of conventional lA, only the 
spatial dimensions of the transmitter and the receiver are used for interference alignment. On the other hand, in the 
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TABLE I 

Required Number of Users to Maintain DoF M (= Nr/2) with R'rf^^ < M 



SNR 


Nr^2 


Nr = 4 


Nr = 6 


5dB 


(3.16)° 


(8.44)" 


(41.3)" 


lOdB 


(10)" 


(8.44 X 10^)" 


(1.31 X 10*^)" 


15dB 


(31.6)" 


(8.44 X 10^)" 


(4.13 X 10^°)" 


20dB 


(100)" 


(8.44 X 10^)" 


(1.31 X lO^'^)" 


25dB 


(316)" 


(8.44 X 10**)" 


(4.13 X lO^'')" 



proposed OIA, additional user dimensions provided by K users are also used to align interference. Therefore, the 
whole spatial dimensions M can be secured and translated into achievable DoF. 

Lemma 4: If we define AR{P) as the difference between the achievable rates for an infinite number of users 
and K users, AR{P) is upper bounded on the average rate loss term (l22l ) as 



AR{P) = R' 



OIA dOIA 



R 



[K] 



^ oloss 



Proof: It is satisfied that 



Elog2 



I. , + ilFH(i)H(i)tFt 



M 



Elog2 



I,, + ZFH(i)H(i)tFt + — y FHWH»tFt 

M ^ M 



M 



i=2 



Elog2 



pa 



M ^ 



i=2 



<Elog2 



P' 



Ij^^ + i_yFHWH«tFt 



i=2 



This is because ^^^2 H^^^H^*)^' is positive semi-definite so that 



i=2 



Therefore, AR{P) is upper bounded by the average rate loss term. 



Theorem 4: To achieve DoF M per user, the number of users per each group, K, should increase as 



V 



i\/3 2<5/^^ - 1 



(25) 
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That is, K should scale like the interference power, P", such that 



K (xP' 



Proof: See Appendix |El 



In Table H the number of users to achieve DoF Af(= ■^) per user is summarized when the average rate loss 
term is set to be less than M (i.e., R^°jf-^ = M) which means a 3dB SNR gap between and R'j'^^ in the high 
SNR regiorQ. When = 2 and a = 0.8, for example, 40 and 101 users are required at 20 and 25dB, respectively, 
to achieve DoF 1 (= Nr/2) per user. For the case that = 4 and a = 0.5, we need 85, 844 and 8438 users for 
15dB, 20dB and 25dB, respectively, to achieve DoF of 2 (= iVj?/2). 

In Fig. 121 the achievable rate per transmitter with OIA schemes is plotted when the number of users is scaled. 
The antenna configuration and the relative path loss are given by {Nt, M, Nr) = (1, 1, 2) and a = 1, respectively. 
With the scaled number of users such as K oc P"*^^, the achievable DoF is maintained as M. This is because 
the scaled number of users achieves a lower rate than an infinite number of users within a constant gap. Now the 
following theorem quantifies how the required number of users can be reduced as the target DoF per user decreases. 

Theorem 5: To achieve DoF M' per user such that 

M - aM < M' < M, 
the number of users per each group, K, should be increased as 




Proof: If DoF M' is achieved using OIA scheme with K users, M' can be represented as 




= M'. 



Thus, we need to find the number of users K satisfying 



= M -M'. 



The rate loss is set such that AR{P) = Rf, 



,OIA 



RlK] <Mlog2{SNR) - Mlog{SNR/2) = M = R\ 



(loss 
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Since AR{P) < Mlogs (l + j^D), 



Mlog2 (1 + tt-D) 
lim ^ „ ^ = M - M', 



P-^oo log2 -P 

and substituting D in above equation yields 



VI. Comparison with Conventional Opportunistic User Selection 

OIA aligns interference by selecting the user whose interfering channels are most aligned. This section compares 
OIA scheme with conventional opportunistic user selection schemes, and the computational complexity of each 
scheme is analyzed. 

A. Maximum SNR User Selection (MAX-SNR) 

In maximum SNR user selection, each user maximizes the achievable rate ignoring the interfering channels. At 
the kth user in the first user group only considers its own channel and finds the postprocessing matrix such that 



Fi,fc = arg min logg 
F' 



r(l)TT{l)t 



(26) 



which becomes Fi^^ = [vi(Ai^fc), . . . , vjv/(Ai^fc)]^ where Ai^^ = ^Hj^ jj . Using this simplification, the corre- 
sponding rate for the kth user is log2 nm=i (l + jf^m{-^i,k))- Thus, the feedback information at the fcth user 
becomes nm=i + ^^m{-^i,k)) ^ and the index of the selected user, ki, becomes 



ki = arg max TT 1 + ir7Am(Ai,fe) . 
i<k'<K ^^^\ M J 



B. Minimum INR User Selection (MIN -INR) 

In minimum INR user selection, each user employs a postprocessing matrix that minimizes the rate loss term, 
C}°^^ in dill. A postprocessing matrix at the kth user in the first user group is given by 

F'/;''^=[vAf+l(Bi,fc),...,V2Af(Bi,fc)lt 



where Bi^^ = h[^],h[^],^ + h[^],h[^],^ Then, corresponding C^°f becomes 

2M 



C}7 = log2 n (l + ^A^(Bi,;i.)) 

m=M+l ^ ^ 
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Thus, the feedback information at user k is nm=M+i + "jf Am(Bi fc)), and the index of user selected at the 
first transmitter becomes 



/ pa \ 

/ci=argmin TT 1 + — Am(Bi^fc) I 
i<k'<K ) 



m=Af+l 

C. Maximum SINK User Selection (MAX-SINR) 

In maximum SINR user selection, the user offering the maximum rate is selected. Each user adopts the post- 
processing matrix to maximize the achievable rate, so the kth user at the first transmitter finds the postprocessing 
matrix to maximize Ci_^: given in (|7]l. 

When M = 1, the postprocessing vector Fi ^ € C^"^^ at the A:th user maximizing SINR such that 

1 1 k = arg max ttt — tttt 

is obtained by solving a generalized eigenvalue problem and becomes 

Fi,fc = I^i}, ■ vi (^(L];^^)1'^Ai,fcL^j 

where Li ^ is Cholesky decomposition of positive semi-definite matrix Itv^^ + -^Bi ^ such that lJ jj^i^k = Itv^ + 
"I^Bi^fc. At this time, corresponding Ci,fc becomes logs (l + '^i ((Ll^fc)^p^Al,fcL^,fc))• 
When M > 2, the postprocessing matrix Fi ^ at the fcth user to maximize Ci^k such that 

i-M+t [jjti^ f^ti^ f^ + AT 2^i=2 ^l,fcWl,fc j * 

Fi k = arg max logs 

F' 



Am + AT l^i=2 * ^l,fc^l,fc * ' 



is hard to find and generally unknown. 



D. Time Division Multiplexing (TDM) 

In this scheme, only one of the three transmitters serves its selected user at any point of time. Therefore, the 
selected user does not receive any interference from other transmitters. Each user finds the postprocessing matrix 
to maximize the achievable rate, so the postprocessing matrix at the transmitter is the same with the postprocessing 
matrix in MAX-SNR scheme given in (l26l ). Also, the feedback information at the user and the user selection 
criterion are exactly the same as those of MAX-SNR scheme. Because only one selected user is exclusively served 
by the TDM approach, the achievable DoF per transmitter becomes 
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TABLE II 

The complexity of various operations for G g C 



Operation 


Complexity (flops) 


qG, G + G 


2mn 


\\G\\f 


4mn 


G®G = GGt 


8mn^ — 2mn 


Gram-Schmidt Orthio. 


Sm-n? — 2mn 


Singular Value Decomp. 


24m^n + 48mn^ + 54n^ 



TABLE III 

The complexity of various schemes 



Scheme 


Complexity 


Ratio 




(A-,JV„^^) 


MAX-SNR 


K X (128iVi -Nl + fiVfl) 


98.4% 


MIN-INR 


K X (ISONl + 3Nl + ^Nr) 


100 % 


OIA 


K X {m'i + 2ND + (i30iv| + ml) 


6.15% 



E. Complexity Analysis 

The computational complexity of each scheme is represented by the number of floating point operations (flops) 
|[20l Chap. 1]. An addition, multiplication, or division of real numbers is counted as one flop, so a complex addition 
and multiplication are counted as two flops and six flops, respectively. For an m x n complex matrix G G c™^" 
(m > n), the flops required for several matrix operations are summarized in Table HI] where the operation Cg) is 
defined as G G = GG^. 

In MAX-SNR user selection, each user requires one (g) operation, a single SVD, 2M real additions and M 
real multiplications to find feedback information. Correspondingly, the total computational complexity becomes 
K X (SNrM^ - 2NrM) + Kx (24Af| + 48iV| + 54A^|) + Kx2,M = Kx (128A^^ - + I^r) flops. In MIN- 
INR user selection, two Cg) operations, two matrix scaling, a single matrix addition, a single SVD, 2M real additions, 
and M real multiplications are required at each user to find the feedback information, so the total computational 
complexity becomes K x 2{8NrM'^ - 2NrM) + Kx 2N'^ + Kx 2N'^ + Kx (24iV| + 48iV| + 54A^|) + Kx3M = 
K X (ISOA'^^ + 'iN'j^ + |iV/j) flops. Note that the postprocessing matrix should be calculated to find feedback 
information both in MAX-SNR and MIN-INR schemes. On the other hand, the proposed OIA scheme requires two 
Gram-Schmidt orthogonalization, two (g) operations, one matrix addition, and a single || • \\p operation to construct 
feedback information. The selected user needs 130A'^^ + 'iN'^ additional complexity to find the postprocessing 
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matrix, Ff^. Therefore, the total complexity of OIA scheme becomes K x A{8NnM'^ - 2NrM) + K x 2N\ + 
K X ml + (130A^| + 3iV2 ) = Kx (8iV| + IN^) + (130iV| + SiV^ ). 

The computational complexity of various schemes is summarized in Table |lll] and plotted according to the number 
of users, K, in Fig. |4]when Nr = 4. The complexity of the proposed OIA scheme is just about 6.15% of MIN-INR 
user selection scheme's complexity when the number of receive antennas, Nr, and the number of users, K, are 
sufficiently large. 

F. Performance Comparison 

In Fig. |5j the proposed OIA scheme is compared with other user selection schemes in terms of achievable rate 
per transmitter when {Nt, M, Nfj) = (1, 1, 2). The relative interfering channel power a, and the number of users, 
K, are set to be 1 and 50, respectively. It is sho wn that the proposed OIA scheme achieves a similar rate to the 
conventional MIN-INR user selection scheme while it requires much less computational complexity than MIN-INR 
user selection. It is also shown that the proposed OIA scheme significantly outperforms the conventional MAX-SNR 
user selection in high SNR region. For a finite number of users, the achievable rates using MAX-SINR, MIN-INR, 
and OIA schemes are saturated in high SNR region. On the other hand, the TDM scheme achieves a DoF of ^ and 
outperforms OIA above 50dB SNR. 

Fig. |6] shows sum rates of the proposed OIA scheme and the conventional user selection schemes for the antenna 
configuration {Nt, M, Nr) = (2,2,4). The number of users, K, and a are set to be 50 and 0.7, respectively. 
Similar trends can be observed to Fig. [5] although the sum rates are greater than those in Fig. |5] due to the increased 
spatial dimensions. As predicted in Theorem |2j the achievable DoF becomes 0.3 (= 1 — a) in high SNR region. 

VII. Conclusion 

We have proposed a novel opportunistic interference alignment (OIA) and analyzed its optimality in terms of 
the achievable DoF in the three-transmitter MIMO IC channels. When Nr = 2M and Nt > M, the proposed 
OIA scheme has been shown to achieve M DoF per user by opportunistically selecting the user whose received 
interference signals are most aligned with each other. Contrary to conventional lA which is known to achieve ^ 
DoF per user in a three-transmitter M x 2M MIMO IC, the proposed OIA scheme does not sacrifice the spatial 
dimensions in aligning interference signals and secures the full spatial DoF by exploiting the user dimensionality 
for interference alignment. Furthermore, the proposed OIA scheme does not require global channel knowledge at 
the transmitters but needs only scalar value feedback from each user for user selection. We have also proved that the 
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full DoF can be achieved even with a finite number of users if the number of users grows with an appropriate scale. 
Finally, the proposed OIA scheme has been shown to have advantages over conventional user selection schemes 
for interference mitigation in terms of both computational complexity and achievable rate. 

Appendix A 
Proof of Lemma [T] 

Using the unitary matrices Y and Z in AA^ + BB^ can be rewritten as 

AA^ + BB^ = AY(AY)^ + BZ(BZ)^ 

M 

= 5] (a,a| + b,b|) . (A.l) 

1=1 

Also, we decompose bj as 

bj = cos OiELi + sin ^jCj, (A.2) 

where 9i is the ith principal angle and e, is a unit vector (||ej|| = 1) orthogonal with such that _L a^. 

From the property of principal vectors, a^ _L a^ and bj _L bj for i / j. Also, from the relationship between the 
principal angle and the principal vector given in ([5]), it is satisfied that 

(AY)tBZ = [ai,...,aM]T[bi,...,bM] =D 

where D = diag{fii, fi2, . . . , /iA/) is diagonal matrix, so aj _L bj for i ^ j. Thus, it is satisfies that span{ai, bj) _L 
span{aij,hj), for i ^ j and equivalently, span{a.i,ei) _L span{aij,ej) for i / j. Using this and from the fact that 
a, _L ej, we can conclude that {ai, . . . , ajvf, ei, . . . , e^/} becomes 2M orthonormal bases of C^*^. 
Using (IA.2I ). b.jb| becomes 

bjb| = (cos OiSLi + sin0jej)(cos OiSn + sin^jej)^ 
= cos^ 6i ■ Rial + sin^ 9i ■ e^ej 
= ^-ajal + (1 - M-)eie|, 

and (lA.ll ) can be rewritten as 

M 

AAt + BBt = (aja| + bjbj) 

1=1 

M 

= J][(l + /.2)aja| + (1-^2)^4 . 

i=l 

Thus, AA^ + BB^ has the eigenvectors {aj} and {ej}, and has ordered eigenvalues 
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Appendix B 
Proof of Lemma [2] 

Consider the unitary matrix Uj^ = H^^^oiaH^^].^ satisfying U|.Ufc = UfcU|. = Iat^. Because the chordal distance 
between two subspace is invariant with a rotation, the chordal distance between two subspaces spanned by two 
interfering channels at the kth user in the first user group satisfies 

'^^(HS,Ha) = ^?(U.HS^l,U,Ha) 
At the selected user kf^^, the chordal distance between two interfering channels can be represented as 

j2^fT(2) tt(3) \ _ ^2/Vt(2) tt(3) n 

72/TT..-tr(2) TT. Tj(3) 



min di(Uk'UYi„\Jk'liW 

k' ' ' 



= min dl{-Hfl^,W). (B.l) 

where Crand = |UiHp], . . . ,UxH[^]^| is a random codebook of size K whose elements are i.i.d. and isotropic 
in Qnr,m{C). After simple manipulations and taking average on both side of dB.ll ). we obtain (TTl . 

Appendix C 
Proof of Lemma [3] 

Because F given in ( [20b minimizes the rate loss term, the average rate loss term i^j^j can be rewritten and upper 



bounded on 



R'(^] = Eh(.) min logs 



pa 3 



M ^ 



(a) 



i=2 
3 



Ia/ + ^ 5] F'H»Ea(.) [A(^)]H(^)tF't 



M 

i=2 

<Eh(., Mlog2 1 + — Yl 

V m=M+l 
(c) / P 

< Mlog2 I 1 + —D 



M 

where C = ^^^2 H^*^H^*^^. In above equations, inequality (a) is because of the independency between H^*^ and 
A^*^ and Jensen's inequality. In the OIA scheme, both the feedback information at each user (Theorem [Hi and 
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the user selection dTSl ) are independent with A^*), so it is still satisfied that E[A(*)] = MI a/. Thus, inequality (6) 
holds by the same reasons with ([TTI)-([T3]). The inequality (c) can be easily shown with Theorem [U ( [T6l ). Jensen's 
inequality, and ( fT9l) . 



Appendix D 
Proof of Theorem [3] 

The term ii^^" in ([D has lower and upper bounds such that 



i?g>Elog2 



I^, + ^FH(i)H(i)tFt 



(D.l) 



and 



i?f;j<Eiog2 



i=2 



Elog2 



+ Z,FH(i)H(i)tFt 



loss 



(D.2) 



Because both FH^^^H^^^^F^ and ^^^2 FH^'^H^'^^F^ are positive semi-definite matrices, the lower bound (ID.ll ) 
is obtained by neglecting the interference term ^^^2 FH^^^H^^^^'F^ The upper bound given in (ID .21 ) is obtained 
because 



I^, + ilFH(i)H(i)tFt + — y FH»H«tF 



M 



i=2 



Im + HWtFt (Im + ^ E FH»H»tFt) 'fH^^) 



i=2 



pa 



i=2 



(a) 
< 



Im + h(i)^f1"FH(^)| • 



Im + — J]]FH(*)H«tFt 



i=2 



Note that the inequality (a) holds since 



HWtFtFHW -H«tFt(lM + ^ J]FH«H»tFt) FH^^ 



1) 



i=2 



is positive semi-definite so that (H(i)tFtFH(i)) is bigger than (^H(i)tFt Im+j^ EL2 FH^^^H^^^tFt 
for all m (1 < m < M). 

In Lemma [3l the rate loss term, -Rj'^p becomes zero as ii' — 00 because limi^_i.oo -0 = 0. Applying the squeeze 
theorem for the two bounds, we can easily show that 



fh(i 



lim <iJ = Elog2 



+ ^FH(i)H(i)tFt 



where F is independent of H^^^ and uniformly distributed in C^^^^" so that the distribution of F is the same with 



the distribution of the postprocessing matrix of a user in the Nr x 4f x M BD system. 
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Appendix E 
Proof of Theorem H] 

As noted in Corollary (TJ the DoF M is achievable as K — )• oo. Thus, if we maintain Ai?(P) within a constant, 
we can achieve DoF M even with a finite K. Using Lemma [3] and Lemma IH we can make Ai?(P) less than 5 
such that 

Ai?(P) = i?f4 - i2p4 

< R'^ 

<Mlog2 + 

< 5. 

Thus, we can achieve DoF M, if the upper bound of rate loss term is maintained within a constant 5 such that 

MIog2 + 

Substitute the main order term of D in dTSl) in above equation, we obtain K to maintain the DoF as in ( |25l) . 
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Fig. 1. The system model considered in this paper is shown above. Each muhiple antenna transmitter selects one user from its group. 




Fig. 2. The achievable rate per user of OIA scheme varying a when {Nt, M, Nr) 



= (2,2,4). 
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Fig. 3. The achievable rate per user of OIA scheme with scaling K (x P when {Nt, M, Nr) = (1, 1, 2) and a — 1. 
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Fig. 4. Complexity of various user selection scheme according to the number of users K when Nr — 4. 
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Fig. 5. The achievable rate per user of various user selection schemes according to the SNR when {Nt , M, Nr) — (1, 1,2), a — 1 and 
K = 50. 
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Fig. 6. The achievable rate per user of various user selection schemes according to the SNR when (A'r, M, Nr) = (2, 2, 4), a = 0.7 and 
K = 50. 



